Study of the Magnetic Film Materials by Horizontal Scanning Mode for the Magnetic 
Force Microscopy in Magnetostatic and ac Regimes. 
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The magnetic force microscopy inverse problem for the case of horizontal scanning of a tip on a 
linear magnetic film is introduced. We show the possibility to recover the magnetic permeability 
of the material from the experimental data by using the Hankel (Fourier-Bessel) transform inverse 
method (HIM). This method is applied to the case of a layered slab film as well. The inverse problem 
related to the ac MFM is introduced. 
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I. INTRODUCTION 

In the theory of inverse magnetic force microscopy 
(MFM), the recovering of the magnetic permeabiUty is 
based on the interaction of a tip with the stray field of a 
magnetic material. So far the inverse problem has been 
solved in magnetostatic considerations for different cases 
such as spherical^, seminfinite^ and slab"^ geometries. In 
order to measure the magnetostatic interaction force, the 
direction of the MFM tip movement has been always as- 
sumed to be perpendicular to the surface of the sample. 

This vertical movement mode of the MFM entails 
a mathematical procedure based on the Laplace trans- 
form inverse method {LTNV^i^i^i^ applied to experimen- 
tal data. Solutions of the LIM are not unique since 
the numerical Laplace transform inversion is intrinsically 
unstable^'^. In other words, small variations in the initial 
conditions or the noise present in the experimental data 
cause large variations of results. 

In this paper, we assume lateral displacements of the 
tip which allow us to consider the use of the Hankel trans- 
form inverse method (HIM)^^, an alternative approach 
to LIM, for recovering the magnetic permeability ^Xr- We 
first present the inverse problem related to the MFM for 
a finite magnetic slab. We extend the proposed Hankel 
inverse algorithm to the case of a finite layered film and 
recover the magnetic properties as well as the geometrical 
characteristics of each layer. 

Finally, we introduce the HIM for the ac MFM regime 
to study a magnetic slab. 



II. SLAB FILM MAGNETOSTATIC INVERSE 
PROBLEM 

Let us consider the tip as a magnetic point dipole m 
with components (0, 0, to^) placed at the position (0, 0, a) 
outside a magnetic slab film of thickness b. The film is 
located at z = parallel to the xy plane. We assume the 
relative magnetic permeability iir — fJ- / fJ-o oi the material 
to be constant, where /xq is the magnetic permeability in 
the vacuum. The scalar potential (f) satisfies the Poisson 



equation 



V^(/) = TO ■ VS{x) Siy) S{z - a). 



(1) 



The general solution of the Poisson equation is given 
by a sum of a particular solution 0i and a homogeneous 
solution 02. This decomposition of (j) corresponds to 
the decomposition of the induction i?(r) (where B{f) = 
— /zoV0), into the sum B{r) = Bi{r)+B2{r), where Bi{r) 
represents the direct magnetic field induction due to the 
dipole, and B2{r) is the refiected magnetic induction. 
For the dipole to, we have 



/io TO • r 



(2) 



where ri = (x, y,z — a). 

Now, we expand 9^1 (r) in terms of the Bessel functions. 
For that we use the following identity 



1 



dfce-'^l^l Jo(fcp)- 



(3) 



Here and throughout, Jg denotes the Bessel function of 



Oth 



order of the first kind and p = x 



y 



Eq.® 



and the recurrence relations of Bessel functions permit 
us to express the z-component of 4>i{r) from Eq.([21) in 
the following^ compact form. 



/>.i = ^/ dke-'^^^—HUkp). (4) 



The distribution of the magnetic field induction B{r) is 
given by 

(Bi{r)+B2{r) z >Q 
B{r) = <^ Bajr) -6< z < (5) 

where Bz{r) and Bi{r) are the penetrating and transmit- 
ted magnetic fields respectively. 
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The following governing elliptic differential equations 
hold for B{r) 



V'B,{r) = 0, (z = l,...,4). 



(6) 



Then, analogous to Eq.Q, we can write the z- 
component of the scalar potentials'^ as follows: 

/•OC 

0.2 - ^ / dfcfc42(fc)e-'=('^+^)Jo(M, 



An 



47r 7o 



^zA = dkkI+{k)e-''^''-'-''^Jo{kp). (7) 

47r Jo 

The coefficients 1^^ (k) can be obtained if we impose the 
continuity boundary conditions on 4>zi and fj, d4>zi / dz at 
the interfaces {z = Q^z = —b). Notice that we do not 
consider the solutions which are divergent as z —^ ±cxd 
for the fields Bzs and Bz4- We are interested in the co- 
efficient Iz2{k) since Bz2 is the only magnetic field that 
interacts with the tip. For this coefficient we have from 
Eq.dll) 



(8) 



where = (e^'''^(^r + l)^^(A*r^l)^)^"^- The expressions 
for the incident, reflected, penetrating and transmitted 
fields respectively are 

Bz, = ^ / dkk'e-'^'^-^^Mkp), (9) 



^ fiornz_ / dkk'lzJpr,k)e-'^^+^^Jo{kpllO) 
47r Jo 



Bz. - 



firuz 
47r 



X e-'=Vo(M, 



(11) 



Bz 



An 



dk e IzAt^r, k) e-'=('^-^-^) Jo(M)12) 



The interaction force between the tip and the slab is 
given by 



Fzia)^-li-lrn.U 



(13) 



and it can be rewritten using Eq. lTOl) as follows 



Fzia) 



47r Jo 



dfcfc^/,2(/ir,fc)e~'"'= Jo(M- (14) 



Let us assume that we need to recover the unknown pr 
from the MFM data (i.e. the force Fz). Consider a hor- 
izontal displacement of the tip such that the distance a 
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FIG. 1: Plot of the simulated force F^i^a^, p) used in this work. 
A noise corrupted data have been added. 



between the tip and the magnetic film remains constant, 
a = ao. 

Recalling the definition of the Hankel transform oper- 
ator 



^^[imip) = / kJoikp)fik)dk, 



we can write Eq. (|M)l as follows 



Fz{ao,p) - ^^J{o[fc'/.2(/ir,A:)e~2''-«] 
An 



(15) 



(16) 



The data inversion process involved in our algorithm re- 
duces to a simpler operation, the so-called Hankel trans- 
form inversion. The forward and inverse transformations 
have the same operational form^° = IKo. 

Let us now describe the algorithm for finding p,r start- 
ing with the experimental force data. The mathematical 
inversion should be performed for the pairs [pi, Fz{ao, pi)] 
where p is the Hankel transformed variable. We apply the 
inverse operator 5£o ^ to both sides of Eq. lTBl) . 

Consider the quantity ^Q^[Fz{ao, p)]{k) to be known 
for some collection of values of the wave-number k. Using 
these data it is possible to determine the pair {pr , b) or 
(Mrj-fz2) as follows. Assuming A: as a parameter, we can 
take any nontrivial pair ki ^ k2 and write the following 
nonlinear system of equations: 



Ti^^^:Ko'[Fziao,p)]ik^)-Iz2{pr,ki) = 
I^^^'K-\Fz{ao,p)]{k2)-Iz2{^lrM) = 0,(17) 



where T(ao, k) = Ane'^'"^" / poml- 

The system (jl7p allows us to obtain the quantity pr 
from the force experimental data. To prove that the 
described method gives the value of we perform 
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FIG. 2: The magnetic permeability /i^ and the coefficient 1^2 
recovered by using the Hankel transform technique. 



a simulation^ii^ with a set of faked MFM data with 
ao = 0.435, b = 0.2, p = 0.2 and fir = 0.1. Here and 
throughout the paper the units are measured in fim. The 
corresponding graph is shown in the Fig.l. 

Fig. 2 displays the results of the recovering of Hr and 
the coefficient Iz2il^r, k). The results agree with the the- 
oretical value of the parameter /i^ used in the forward 
problem shown in Fig.l, in a large range of values of the 
wave-number k. 



III. LAYERED FILM MAGNETOSTATIC 
INVERSE PROBLEM 



The inverse problem for the slab can be generalized in 
order to study the magnetic multilayered films. 

Let us consider a finite slab film that geometrically 
can be viewed as a collection of M layers. Denote the z- 
coordinate of the j-th layer by zj . We say that the lowest 
layer is designated as layer 1, while layer M is the top 
of the slab and zm = 0. The slab has a total thickness 
zi = —b. We assume a constant magnetic permeabilities 
l^rj {j — 1, M) to be different for each layer. 

The distribution of the magnetic field B{r) can be writ- 
ten as 

Biir} + B2ir) z > 0, 

B3{r) Zj < z < Zj+i, (18) 
i?4(r) z < zi. 

The elliptic differential equations V^Bi = (i = 
1, . . . , 4) hold. The z-components of B2, B3 and B4 are 



given by 

Bz,, = 



X e 



Bz, = 



lk)mz 

47r 
fJ.mz 

-^•Vo(M, 

fiom 



dfcfc2/(M+i)(fc)e-'=('^+^)jo(fcp)j(19) 



dke[Ii+ik)e'^^+p-ik)e 



— kz 







(20) 



47r 



dfcfc2j^,(fc)e-'=('^"^-^i)jo(M-(21) 







The coefficients Iz^ (k) may be obtained by imposing 
continuity boundary conditions at the planar interfaces 
Zm = and z — Zj on Bz^ in a way similar to that 
described in section II. 

Taking into account Egs. p^ and p^ . we can write 
the expression for the force between the tip and the lay- 
ered film as follows 



47r 



dkk'li^'+'^e''^' 



Mkp). (22) 



Notice that Iz2^^^^ contains all the geometrical and 
physical information related to the layered film. This 
coefficient is a function of the permeability and position 

of each layer, i.e. li2^~^^^ — ll^'^^\k, Zj, prj)- 

Taking into account the definition of the Hankel 
transform Jfo, we apply its inverse to both sides of 
Eq. ((22|) . For M couples of values kj+i ^ kj we get the 
following nonlinear system, which is analogous to p7|) : 



T{ao,k 



M+l, 



k^ 



J{o'[Fz{ao,p)]ikM+i) 



^1 



-J{o'[F,(ao,p)](fcj 



iir'\k,,^,) = 

... (23) 
-lT'\k,) ^ 0. 



In order to get the force noise corrupted data, we have 
followed the simulation procedure used in section II for 
the slab geometry. Fig. 3 shows the result of applying the 
described inverse method to a film having three layers. 

Notice that the magnetic layered film-tip interac- 
tion problem involves an attenuating magnetic field 
penetration^. If the experimental data of the force reflect 
this fact, then better results in the inverse procedure for 
layers situated next to the top of the slab {z = 0) can be 
archived. 



IV. ac REGIME FOR A SLAB FILM 

It is well-known that one achieves better results using 
the MFM if one considers the ac regime rather than a 
magnetostatic interaction between the tip and the filin^. 
The dc measurement is easily corrupted by noise, such 
as vibrations and 1// electronic noise. Much higher 
sensitivity of MFM can be achieved in the ac mode by 
driving the cantilever at its resonant frequency by that 
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FIG. 3: Recovery of the magnetic permeabilities /iri, Mr2 and 
fir3 for a slab with M=3. The theoretical values assumed for 
firj are 0.41,0.9,0.71, respectively. 



Replacing a by a' in ([9]), we get 



Bzi — 



47r 



47r Jo 

/•OO 

Jo 

(24) 

The following similar expressions can be obtained for Bzi 
= 2,3,4): 



Bz2 = 



47r 



X Jo(fcp) + ^^^'2 cos2(7«dt) / dfcfc4e-'=(''+^) 



X Iz2ik,fi')Joikp), 



(25) 



one reduces environmental noise and increases the sig- 
nal to noise ratio using lock-in or frequency modula- 
tion techniques^. To develop and to solve the inverse 
problem in the case of the ac regime with the inclusion 
of losses we follow the method developed by Badia in 
Ref.Q. 

First, let us assume the tip (0, 0, tUz) to be harmoni- 
cally driven at a constant height above a magnetic slab. 
The cantilever oscillation amplitude is considered to be 
constant as well. This condition is used for dissipation 
imaging, an imaging technique which is sensitive to non- 
conservative interact ion a ^ ^ i ^ ^ . 

We consider a small vertical oscillation of the tip 
around the point located at the distance a from the film. 
Denote by a' the distance between the tip and the film. 
Then we have the relation a' — a + A' cos{wdt), where 
A' is the amplitude of the tip oscillations, which satisfies 
the condition A' <^ a. The quantity Wd is used for the 
driving angular frequency of the MFM tip. 

Taking into account the nondispersive ohmic relation 
J = anE and the well-known time-dependent Maxwell 
equations, we get the following relation for the magnetic 
field inside the film V^B = pCndB / dt. The quantity 
(T„ is the conductivity and p is the permeability of the 
material. 



The relative permeability /i^ to be determined is re- 
lated to the real part fx' of /2 as follows: ^' = ^rMOj where 
fi — {p! + ijJ'"}- Now we are going to determine /i as a 
first step in order to get /z^. To that end, we shall find 
the force acting on the tip and solve the inverse problem. 



/•OO 

Bz3 - dkk'[I+{k,n')e''^'+i;,{k,n')e''^]e~'^ 

poo 

X Ja{kp)-—A'cos{wdt) / dkk^[I+{k,p')e^' 
471- 7o 

+ /73(fc,M')e'le-'=Vo(M + "^A'^ cos" {wdt) 



OO 



X / dkk'[I+{k,p')e'^-+i;-,{k,p')e ^ 
lo 

X e-'=Vo(M, (26) 

Bz, - ^ / dkk'e-'^'^-^-'>hz,ik,Mikp) 
in In 



A'coa{wdt) / dkk^^e'^^^-'-^hzAik.p') 



in 

X JM + cos^ (wdt) r dfcfc4e-'=('^-^-^) 



X Iziik, ^i')Joikp). 



(27) 



Notice that the coefficients Izi{k, fj,'){see Eq.®) are sim- 
ilar to the ones found in the magnetostatic caseS.. 

In order to find the coefficients /^(/c,^'), we con- 
sider the following representation of the incident, re- 
flected, penetrating and transmitted fields as Bzi{t) = 
Bit + K^Wd, t) + Bl^^iwd, t), where (i = 1, . . . , 4). 

Let us denote by w a formal variable which takes 
the values Wd and 2'Wd- Assuming the notation 
i3'"^=Re[i?e'"^] and taking into account the time- 
independent solutions p fT^ for B^f (i=2,3,4), we get 
the following governing differential equations 

v'B:,{t) = 0, 

V'B:,{t) = i25-'\w)B:S), (28) 

y'B:,{t) = 0, 



5 



where the skin depth S{w) is defined by S{w) = 
{2/wfla„y/\ 

We note that Eq. ([28|) for is analogous the equation 
V^B = (1/A^)_B, which arises from the London equa- 
tion for the superconductor (A is the London penetration 
depth). 

In a similar way to Eqs.([9l fT2|) . we find solutions B^f 
and S^f"^ {i = 2,3,4) of system 1^ as well as B^f in 
the form: 



B 



47r 



+ 



A'^ / dkk^e-''^''+''h,2{k,ti').h{kp), 



TjWd _ 

J=>z2 - - 



47r 



(29) 



A' / rffcfc3e-'=('^+^)27/(fc,^')^o(M, 

(30) 
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FIG. 4: The graph shows the real part of the magnetic per- 
meability Hr and the thickness of the film b by HIM applied 
to ac MFM. 



and 



and 



= iS^^'Vo dkk^e-''^'^+^'>ll-^{k,f,')Mkp) 



The corresponding expressions for the penetrated fields 
are: 



B 



dc 
z3 



47r 



dkk'[I+(k,p')e''^+I~(k,p')e' 



X e-'^'Mkp) 

poo 

+ ^A'^ dkk^[lUk,p')e'^^+I-,{k,t/)e-^^] 
X e-'^'Vo(M, (32) 



and 



Xr3''"(^-M')e-'le-''-Jo(M, (33) 

poo 

-I- ilr'{k.i^')e~'^']e~'^-Mkp). (34) 
and the transmitted fields are given by 



(31) 4T(fc,/^') = -^A'J dfcfc2e-'^(— ^)I,7(^,A^') 
X Jo{kp), (36) 

/•OO 

4'r(fc,M') = / dA:ft^e-^-(«-^-'')l,^-^(fc,M') 

X Jo(M- (37) 



To obtain the coefficients X^''(fc,^') and X^^''(fc,^'), 
we must impose the continuity boundary conditions on 
each field components at the interfaces z = and z — ~b. 

Now we shall derive the time-dependent force acting 
on the oscillating tip. In fact, in the dipole limit, the 
instantaneous value of this force may be calculated from 
the expression Tz{t) — —da'{ — {l/2)m ■ i?2[a'(t)]}. 

The magnetic force J^z{t) can be decomposed into the 
sum Tz{t) = Ff"^ + F^'^{t). The linear approximation 
with respect to A' for the components of J^z{t) can be 
written as follows 



Bii{k,p') 



47r 



dk k^e-''^''-''-''hzi{k, p)Jo{kp) 



+ / dkk^e-^'^--'-'^h,^{k,^l') 

167r ./n 



X Jo(M> 



(35) 



F. 



dc 



8^ JO 



dkee-^^''Iz2{k,p,')Mkp) 



+ dkk'e-''^^\i:,^\{k,p') 

X Jo(Mcos^"'^ (38) 
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and 



F" 



X A' cos{wdt) 



dk k^e-^'"'I,2ik,fi')Joikp) 



An 



dkk^e-^'"' 



X A' cos{wdt + ifi'^") 



dkk^e-^''" 



|ir2l(fc,M')^o(M 



\j'--\ik,f,')Mkp) 



87T JO 

X A' cos{2wdt + Lp^'"''), 
were 2,7 = ll^^'le"^^" and X.'a"" = &|e'^' 



(39) 



and (p^ 

is the phase lag between the tip-sample force and the 
displacement of the tip. 

In the lowest-order approximation to the ac problem, 
we can assume that is a small parameter—. 

Taking into account that the force is recovered in the 
form F^'= = RelF^e'"""*], we can rewrite the Eq.® as 
follows 

p^oA'tnl 



F, = 



27r 



dkk^e-^''''I,2{k,fi',Wd)Jo{kp), 



where we define Iz2{k, ij! ,Wd) 



{l/A)h2{k,^l') 



(40) 



(l/2)X,7(fc,M',i«rf) + (l/4)j2»^(fc,/.',2u;,). 

Now, we shall recover /i^ and h from MFM measure- 
ments in ac modes. First, recall that the complex am- 
plitude of the ac force F^ = Re[Fz] + ihn[Fz] may be 
experimentally found. 

Here, we show that a straightforward extension of the 
vector inversion method to complex variables can be 
done. Eq. (j40|) can be written using the operator JCq as 
follows 

F,(ao,p) = -^^4^^o[^'e^''""^-2(fc,M',^«d)]- (41) 

Then, we can perform the mathematical inversion for the 
pairs [pi, Fz{af)^ Pi)]. In particular, one can apply the 
inverse operator J{q" to both sides of Eq. (|4T|) and obtain 
the following system for any fci ^ k2 

g{ki) 



kl 

Q{k2) 



1,3 
'^2 



■Ko\h{ao.p)]{k^) 



%o\Fz{ao,p)]{k2) 



Iz2{f^',b;ki) = 
Iz2ifJ-',b;k2) = 0, (42) 



where g{k) = e^'"'° 2Ti/pQA'ml. 

To solve the system for the quantities {p,r,b), 

we use the vector inversion procedure together with the 
MuUer's method^^. The MuUer's method can be used to 
find zeros of a function and can be applied to complex 
value functions''^. This method is a generalization of the 
secant method in the sense that it does not require the 
derivative of a function. 

Fig. 4 displays the results of recovering the magnetic 
permeability and the thickness of the film. These quan- 
tities were considered constants in the process of simula- 
tion using noise corrupted fake data. 




Lateral distance p 



FIG. 5: Plot of the simulated frequency shift A/ versus p. A 
noise corrupted data have been added. 



A. Frequency oscillation inverse method 

In this section we describe the method of solving the 
inverse problem in terms of the shift frequency. 

It is convenient to develop the inverse method using the 
observed frequency shift A/ of the oscillating cantilever 
rather than the magnetic force, acting on the tip. If 
the tip oscillation amplitude A' is small compared to the 
distance between the tip and the slab aq, the relation 
between the gradient of the force dFz{z)/dz and A/ is 
given hy^ 



Af _ 1 dFz(z) 
fo 2k s dz 



(43) 



where /o is the unperturbed resonance frequency and kg 
is the spring constant of the force sensor. 

We can rewrite the expression for the force given by 
Eq. (gO]) in the form 



Fz{z) = r dkk^e-''(^+'^^Iz2{k,p',Wd) 



X Jo(M- 



(44) 



In the case of the horizontal displacement of the 
z— direction oscillating tip, Eq. (|43p can be written as 



A/ = ^""^^f" dkk'e-'^^^+'^hz2{k,n',Wd)Mkp). 

(45) 



Finally, applying the operators and IKq ^ to the 
Eq. (|45|) and considering z — a and k as parameters, we 
get the following system of equations: 



poA'ml 



:Ko^[Afia)]{h) ~Iz2ih, p',wd) = 0, (46) 
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Thickess of the film b 



kg/k^ 

FIG. 6: Plot of an example simulated force F^i^a^, p) versus p 
used in this work. A noise corrupted data have been added. 



bols). The values of a = 0.8,6 = 0.1, /i^ = 0.01 and 
Wd = IQ^rad/s were used for our faked data. 

Fig. 6 displays results of recovering the magnetic per- 
meability and the thickness of the film. These quantities 
were considered constant in the process of simulation us- 
ing noise corrupted data. 



V. CONCLUSION 

Throughout this work, we have developed a theoretical 
background for the recovery of the magnetic permeabil- 
ity in magnetic films by MFM experiments. In addition, 
it has been shown that h may be recovered from the mea- 
surements of the vertical force in the case of a horizontal 
scanning of the MFM tip. 

Finally, we have suggested a solution of the inverse 
problem using the shift frequency of the system. 



where i = 1, 2. 

We have arrived to a nonlinear system analogous to 
(|42p , but this time it is written in terms of the frequency 
shift A/. The system (pS)) can be solved in the same way 
as system p2|) . 

Fig. 5 shows the theoretical shift frequency A/ (con- 
tinuous line) as well as a simulated measurement (sym- 
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